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[B] Force resultnat-moment coupling matrix (Eqn. 3.16) 

[C] Material elastic constant tensor (Eqn. 3.1) 

[D] Laminate moment stiffness matrix (Eqn. 3.15) 

d Diameter of cutout (1.00 inch) 

E Modulus of Elasticity (tension and compression) (psi) 

E Strain gage excitation level (volts) (Eqn. 4.1) 

|f| Force vector (lbf) (Eqn. 3.29) 

h Thickness of a laminate ply (inch) 

G Shear modulus of elasticity (psi) 

K Stress concentration factor 

[K] Finite element stiffness matrix (Eqn. 3.29) 

k Radius of curvature (Eqn. 3.17) 

1 Panel length (in) 

|m| Moment vector (Eqn. 3.13) 

Stress resultant vector (Eqn. 3.12) 

PD Power density (watt/in 2 ) (Eqn. 4.1) 

R Resistnace (ohms) (Eqn. 4.1) 

r Radial distance from the origin (inch) 

[Q] Reduced laminate stiffness matrix (Eqn. 3.5) 

[Q‘] Reduced transformed laminate stiffness matrix (Eqn. 3.10) 

S Applied stress (psi) 



17 



S Compliance (Eqn. 3.23) 

[T] Rotation transom matrix (Eqn. 3.9) 

w Panel width (in) 

X,Y,Z Rectangular 

{5 f Displacement vector (Eqn. 3.27) 

Direct stress 
< Direct strain 

a Applied far-field normal stress 

v Poisson's ratio 

£ d/[2(d/2 + a. ) ] 

lie microstrain (*10^) 

7 Shear strain 

Subscripts 

avg Average 

c Compression 

i,j,k Indices of summation 

1 Lateral direction (parallel to load line) 

max Maximum 

n Notched panel 

a Stress 

e Strain 

t Tension 

t Transverse direction (prependicular to load line) 

u Unnotched panel 

ult Ultimate strength (indicating total failure) 

1,2 Directions parallel and perpendicular to principal fiber 

direction respectively 

oo Infinite panel width 
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Abbreviations 



ASFC 

DOF 

Eps-X 

Eps-Y 

Eps-XY 

FAPF 

FEA 

G/Ep 

Ksi 

LEFEA 

LEFM 

rasi 

NDI 

SCF 



Average Stress Failure Criterion 

Degrees of freedom (at a node) 

Strain in the x direction (€ x ) 

Strain in the y direction (6y) 

Shear Strain (£._.) 

A Y 

First Audible Ply Failure 
Finite Element Analysis 
Graphite /Epoxy 

Thousand pound force per square inch 
.Linear Elastic Finite Element Analysis 
Linear Elastic Fracture Mechanics 
Million pound force per square inch 
Non-Distructive Inspection 
Stress Concentration Factor 
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I . INTRODUCTION 



The ratio of strength to weight is one of the principal means 
of determining the efficiency of a structure. In aerospace appli- 
cations this comparison can be the most meaningful measure. The 
dilemma of designing an airframe for both strength and lightness 
has been with us since the days of daVinci. The quest for ever 
higher ratios of strength-to-weight has led to the development and 
use of high modulus advanced composite materials, principally 
graphite or carbon fibers bonded together in a polymer matrix, in 
the place of metal. 

Major airframe structural components such as wings or bulk- 
heads require cutouts for bolted or riveted attachment, access to 
interior space and passage of control and fuel lines. Timoshenko 
and Goodier [Ref. 1: pp. 78-84], among others, point out that such 
holes in load-bearing structures act to greatly increase the local 
stress and to reduce ultimate strength. This characteristic is 
referred to as the stress (or strain) concentration factor (SCF or 
K). It seems the SCF may have several definitions, depending on 
the material and the researcher. In this report it shall be 
defined as the highest plane strain existing around a cutout 
divided by tne far-field strain; generally called the gross SCF or 
Kgross* Taking into consideration Saint-Venant ' s principle, the 
far-field strain is assumed to be equal to the strain which would 
exist in an ideal, thin, stressed infinite plate if a cutout was 
not present. Stress and strain concentrations, while inextricably 
linked in elastic materials, are not the same. However, since in 
the application discussed here, there is little appreciable 



21 



numerical difference, the term "SCF" will be used to indicate 
either the stress or strain concentration factor. 

When holes or cutouts are necessary in a structural component, 
airframe designers generally have the choice of accepting either a 
significantly lower ultimate load or greatly increasing the compo- 
nent's strength, and thus its weight. In either case the ratio of 
strength-to-weight is reduced in proportion to the highest SCF 
existing within the member. Properly designed ductile metal 
structures mitigate the effects of SCF by plastically deforming 
under high load conditions. This response delays ultimate 
failure, but can also lead to unacceptable reductions in both 
stiffness and fatigue life. 

The metals used in aircraft construction, principally aluminum 
and titanium, can almost always be considered isotropic (many 
manufacturing processes, however, introduce some minor directional 
properties). The magnitude of the orthogonal strains (X, Y and 
shear) existing at a point in a plane isotropic panel result from 
the orthogonal stress resultant at that point and are in propor- 
tion determined by the elastic modulus and Poisson's ratio. An 
applied in-plane stress on an isotropic plate will not induce 
curvature other than, of course, the possibility of the plate 
buckling under compression. Composite plates are termed "quasi- 
isotropic" when they are composed of anisotropic or orthotropic 
laminae stacked with the directional properties arranged in a 
manner to react identically to a true isotropic material to both 
moments and inplane loads. 

Composite laminates typically lack the ductility of metals. 

The high-modulus graphite/epoxy (G/Ep) fibers in general use in 
the aerospace industry allow approximately 1% strain (10,000jU ) 
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in tension and compression to complete failure. Depending on the 
fiber orientation, panels constructed of laminated advanced 
composites with notches or cutouts can demonstrate from sligntly 
less to much more sensitivity to holes or cutouts than otherwise 
identical isotropic metal panels. As shown by Rybicki and Hopper 
[Ref. 2: pp. 15-27], among others, this sensitivity principally 
depends on the type of weave and the orientation of the plies in 
the laminate; that is, it depends on the degree of orthotropy. 

The inherently brittle nature of advanced composite materials, 
their characteristically low strain to failure, coupled with 
manufacturing limitations make their design a far more demanding 
task than that for metals. Other characteristics, however, 
including fatigue and corrosion resistance, light weight, and 
easily tailored directional properties make the design of 
composite structures very attractive, particularly to the 
aerospace designer. 

A- OBJECTIVES AND SCOPE 

This study was designed to investigate the effect of relative- 
ly simple co-cured reinforcement of a cutout on the strain field 
and failure behavior in G/Ep honeycomb sandwich panels subjected 
to uniaxial compression. Honeycomb construction allows very light 
yet exceptionally stiff structures. The objective was to 
determine if a simple and inexpensive reinforcement geometry using 
small volumes of co-cured G/Ep lamina near the cutout could 
significantly reduce local stress concentrations and increase the 
ultimate failure strength in the honeycomb laminate. The idea of 
local reinforcement around holes is not new; Timoshenko noted 
[Ref. 1: p. 82] tnat "reinforcing rings" could decrease the SCF in 
plates with cutouts. The point was to examine the reaction of an 
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advanced composite, a material whose characteristics differ 
markedly from those Timoshenko addressed. 

The research was undertaken with the manufacturer principally 
in mind. Complex or exceptionally thick reinforcement geometries 
are difficult and expensive to manufacture cost-effectively or 
with a high degree of quality assurance. This research used only 
very thin (maximum thickness: 0.028") ply reinforcement in three 
relatively simple geometries. Since facesheets with reinforcing 
plies on both sides would require machining a precise shallow 
depression in both the face of the the honeycomb core and the 
surface of the layup plate, each difficult and expensive tasxs, 
reinforcement was restricted to the outside surface of each 
facesheet . 

This study was limited to one panel size (10.00" x 8.50"), a 
single loading condition (uniaxial compression) and three rela- 
tively simple reinforcement geometries. The 1.00 inch diameter 
circular cutout was reinforced with concentric co-cured round and 
square G/Ep plies around the hole and stiffening strips displaced 
0.50 inch (1 hole radius) laterally from the cutout edge. The 
total amount of reinforcement used varied from 1 to 5 times the 
G/Ep removed from the cutout. Reinforcement was either one or two 
plies symmetrically applied to the outside of both facesheets of 
the panel. A honeycomb core was used, as it would be in an actual 
application, to increase the panel bending stiffness and thus 
eliminate the buckling of the whole panel as a mode of failure. 

The basic panel facesheets were four layer [0,+45,90] 
HMF330C/34 G/Ep cured to a thickness of 0.056 inch. Cured sheets 
were bonded to both sides of a 0.50 inch thick fiberglass/phenolic 
honeycomb core using 3 M , Inc.'s AF-126 (250°F) cured adhesive. 
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The result was a very light, thin quasi-isotropic laminate: 



[0,+45, 90, core ] s with a great resistance to bending. The 
HMF330C/34 is a woven, high- temperature epoxy (350°F) G/Ep fabric 
manufactured by the Fiberite Corporation of Winona, MN. In order 
to reduce the number of design variables the principal axis of the 
reinforcement plies was oriented only in the direction of the 
applied compressive load. This theoretically gives tne highest 
stress concentration and could be considered the worst case. 

B. REVIEW OF LITERATURE 

1. Background and Historical Research 

The subject of notch-induced stress concentrations in 
plates has been extensively documented. The effects reinforcement 
have on the SCF in plates have received considerably less atten- 
tion. Early research concentrated on metals (isotropic materials) 
and focused on defining the stress and strain fields around 
circular and elliptic cutouts. Recent research has been primarily 
in characterizing the response of orthotropic and anisotropic 
materials . 

Kirsch [Ref. 3] is commonly cited as the first to 
determine exactly the stress concentration factor of a cutout in 
an isotropic material from the theory of elasticity. Howland 
[Ref. 4] applied the solution to Airy's equation in polar 
coordinates to determine the magnitude cd: the SCF. One of the 
earliest papers addressing reinforced holes was by Levy, Woolley, 
and Kroll [Ref. 5]. They investigated the effect of both 
reinforced and unreinforced holes on the buckling strength of 
square isotropic plates. They determined that presence of a hole 
caused only a relatively minor reduction in the buckling 
(ultimate) load. 
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A thorough theoretical, closed-form mathematical treatment 
of anisotropic materials with stress concentrations can be found 
in the work of two Russian applied mathematicians, S.G. Lekhnit- 
skii and G.N. Savin. Lekhnitskii [Ref. 6] principally addressed 
the distribution of stress around the edge of variously shaped 
cutouts in unreinforced anisotropic plates and shells under a 
variety of loading conditions. He determined that a plate with 
high anisotropy, as found in strictly unidirectional fiber con- 
struction, could produce a stress concentration factor near 9 when 
the load was parallel (0°) to the principal fiber direction, and 
slightly more than 2 with the load perpendicular (90°) to it. It 
must be pointed out that, in composite materials, the SCF may not 
have a exactly proportional effect on the reduction in the ulti- 
mate strength of the plate. Due to the composite's ability to 
redirect the load path once fibers are broken or lose stiffness 
through matrix degradation, the ultimate strength is not degraded 
as much as would be expected by the presence of the stress concen- 
tration. This phenomenon is discussed in more detail in section 
II. D. 4. 

Savin [Ref. 7] treated the stress and strain fields in a 
plate resulting from a cutout. He addressed the SCF as a function 
of a plate's linear material properties, ply orientation and 
stacking, and its loading. Hole size, reinforcement and geometry 
were not addressed. A computer program was developed by Garbo and 
Ogonowski of McDonnell-Douglas [Ref. 8: Vol. 3] which computes the 
stress and strain field around a cutout based on Lekhnitskii ‘ s and 
Savin's analyses. It was modified by the author for the IBM 370 
and is listed in Appendix P. 
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Substantial research in stress concentrations in composite 
plates was done by Greszczuk [Ref. 9]. He developed a theoretical 
solution for failure stress and stress concentrations in both 
orthotropic and anisotropic material under tension. His method 
was based on the Hencxy-Von Mises distortion energy method, and 
gave both magnitude and locations of the ultimate stress. Rybicki 
and Schmueser [Ref. 10] investigated the effects of laminate 
stacking sequence, lay-up angles, fabrication temperature and 
thickness on panel stress concentrations using finite element 
analysis . 

There is relatively little research into the effects of 
reinforcement around holes in composite plates which has been 
reported in the open literature. Virtually nothing is available 
on the behavior of notched reinforced plates in compression or on 
the effect on the type of failure of using honeycomb in such 
structures . 

Kocher and Cross [Ref. 11] demonstrated experimentally 
that titanium, graphite and steel reinforcement around a circular 
cutout in a composite plate could reduce the SCF and increase the 
ultimate failure load in tension. Their results, however, were 
based on relatively complex, thick reinforcement geometries that 
have not found acceptance in aeronautical design. 

A novel cutout reinforcement method using oonded hoop- 
wound G/Ep disks was addressed by McKenzie [Ref. 12]. The disks 
were used to reinforce both aluminum and G/Ep plates under tensile 
loads. The method proved effective both in reducing the stress 
around the cutout and increasing the plates' strength. 

The team of Daniel, Rowlands and Whiteside [Refs. 13-17] 
did extensive experimental work in characterizing the effects of 



27 



cutouts on a variety of composite materials. They did some 
limited testing of reinforced specimens in tension and found 
proportionately reduced SCFs and increased strength. They 
determined that interlaminar deformations occurred in the boundary 
region of the cutout, an area they defined as extending about one 
laminate thickness from the free edge. This deformation was 
very nonlinear and could cause delamination at relatively low 
loads. Strain levels next to the cutout, prior to failure, were 
found to be higher than the ultimate failure strain of panels 
without cutouts. Based on that, they determined that the SCF did 
not necessarily produce a proportional reduction in strength. 

They recommended keeping the reinforcement close to the hole, 
using stepped diameter plies ("wedding cake") to facilitate the 
load transfer, and using 45° plies in the reinforcement where 
possible. 

Knauss, Starnes, Ilenneke [Ref. 13] tested unreinforced 
0.15 and 0.24 inch thick T300/5208 panels in compression for 
unbuckled and postbuck led strength. They found that under hign, 
but less than normally ultimate stress levels, the laminate around 
the hole could buckle locally, delaminate and initiate total panel 
failure. A micro-mechanical failure mode was postulated where 
limited finer buckling at the point of stress concentration caused 
by local imperfections such as voids, matrix cracking or poor 
fiber-matrix bonds led to total failure. 

2. Summary of Recent Related Research 

This report is the fourth in a series of investigations on 
the character of stress concentrations in composite plates made of 
HMF330C/34 G/Ep. This particular material was chosen because of 
its current use in both the Trident submarine launched ballistic 
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missile (SLBM) and the prototype Lear Fan propjet aircraft and the 
fact that it has a relatively small data base compared to other 
G/Ep prepreg material currently in aerospace use. This research 
was funded by the Strategic Systems Project Office of the Naval 
Sea Systems Command and greatly assisted by Lockheed Missiles and 
Space Co. (LMSC), Sunnyvale, CA. 

The initial project was undertaken by Herman [Ref. 19] who 
investigated the pre- and postbuckled strength of HMF330C/34 
panels loaded strictly in shear. He used a molded-in 45° flange 
around the cutout to add strength to the shear web. He determined 
that this reinforcement method was well-suited to adding stiffness 
to panels that were not buckled but that the panels did not see a 
significant increase in ultimate strength once buckling had 
occurred. 

O'Neill [Ref. 20] demonstrated that reinforcement of only 
one face of a notched panel under tensile loading provided limited 
additional strength. Initially, the reinforcement of only one 
side of a cutout was considered attractive since only a small 
additional manufacturing effort was required. Asymmetric rein- 
forcement, however, displaces the midplane of the laminate (under 
the reinforced area) toward the reinforced side. Uniaxial tension 
tends to pull this local midplane toward the load line, causing 
out-of-plane bending at the hole, wnich results in high shear 
stress between plies, delamination and premature failure. The 
delamination counteracts most of the decrease in stress concentra- 
tion provided by the reinforcement. 

Pickett and Sullivan [Ref. 21] and Bank, et al. [Ref. 22] 
continued O'Neill's research examining tension panels with 
symmetric reinforcement. They showed that suitably designed 
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reinforcement which was symmetric along the axis extending through 
the panel's thickness could both reduce strain concentrations and 
proportionately increase the ultimate strengtn. No delamination 
was noted in their test panels. 

The work reported here extends the idea of symmetric hole 
reinforcement to compression specimens with a honeycomb core. 
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II. APPROACH TO THE PROBLEM 



A thorough investigation into the effect of reinforcement 
around stress concentrations in composite plates must examine 
various materials, hole sizes, panel and reinforcement layups and 
geometries as well as the means and directions of load applica- 
tion. This research addressed only a small portion of the total 
problem. The material, hole and panel size, layup and loading 
method remained constant — only the amount and the shape of the 
reinforcement was varied. Reducing the number of design variaoles 
to only two allowed an analysis of the sensitivity of the SCF to 
certain thin reinforcement geometries. 

A. METHOD OF INVESTIGATION 

To investigate the effects of co-cured reinforcement around 
cutouts, linear elastic finite element analysis (LEFEA) was 
employed to determine the strain field in each panel configura- 
tion. Plots were drawn of the strains existing on a line from the 
point of highest stress concentration at the cutout across the 
middle of the panel (the X axis) and around the cutout and con- 
tours of the three inplane strain fields (Y, X, and shear). These 
are included as figures in the appendices for each geometry. 
Specimens of each configuration were manufactured, instrumented 
with strain gages and finally loaded in compression to failure. 

The analytical and experimental results were compared and the 
failure mode of each panel evaluated. 

Total failure, in this report, is assumed to be facesheet 
delamination, separation and buckling with massive fiber failure 
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such that the panel could not withstand a full reversal of the 
load. Partial failure was facesheet delamination and separation 
without the massive fiber failure and infers that there could be 
significant tensile strength remaining. 

3. COORDINATE SYSTEM 

There are several different right-hand coordinate systems that 
have been used in the analysis of laminated materials. Analysis 
is, of course, independent of the the system used, but more than 
one student has lost his way attempting to compare methods or 
results expressed in different systems by various recognized 
authorities in the field. 

Tne data presented in this report is based in a cartesian 
system with the plane of the laminated panel aligned in the X-Y 
plane. Individual ply orientations are considered to be rotated 
counter-clockwise from the X axis an angle of theta ( Q ) degrees. 
These plies in the layup are assigned a local orthogonal coordi- 
nate system designated the 1-2 axes. The 1 axis, also referred to 
as the principal axis, is considered to be in the fiber direction 
with the highest elastic modulus (E-^). 

Figure 2.1 shows the upper right quadrant of a typical panel 
in the X-Y (global) coordinate system as well as the ply 1-2 
(local) coordinates. This coordinate system was used by R.M. 

Jones in Mechanics of Composite Materials and Ashton, Halpin, and 
Petit [Ref. 23]. Tsai and Hahn [Ref. 24] chose instead to fix the 
X-Y axes to the ply and the 1-2 axes to the panel. The principal 
researchers in the field do not use the same system. 

The panel is oriented so that the area of greatest interest, a 
horizontal plane bisecting the circular cutout, is aligned with 
the X axis, where y = 0.0". The origin is assigned to the center 
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Panel Coordinate System (X-Y) 




Ply Coordinate System (1-2) 

Figure 2.1 Panel and Ply Coordinate Systems. 
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of the circular cutout. The compressive load is applied to the 
panel 90° to this plane, parallel to the Y axis, and referred to 
as O^* The ^ axis is centered at the midplane of the panel and 
extends through the thickness toward the viewer, completing a 
right-hand coordinate system. 

C. SELECTION OF TEST SPECIMEN CONFIGURATION 

The dimensions of the test specimens were chosen to approxi- 
mate, at least in order of magnitude, a typical honeycomb panel 
with a cutout found in many aerospace applications. The overall 
size was limited by the size of the test machine and compression 
test frame. 

Hong and Crews [Ref. 25], among others, demonstrated that the 
stress concentration in orthotropic composites under uniaxial 
loading was dependent on the ratio of hole diameter to panel width 
(d/w). Whitney and Nuismer [Refs. 26 & 27] pointed out that the 
absolute hole size had a significant effect on the stress gradient 
and ultimate strength when the hole diameter (d) was less than 1.0 
inch. 

The cutout's 1.00 inch diameter was chosen, therefore, to 
limit, as much as possible, hole-size effects. The panel was then 
designed as large as practical to reduce the effect of finite 
panel dimensions and still fit into the test frame and machine. 
Hole-size and finite-width effects are addressed in more detail 
in Section II. D. 3. The specimen size, 10.00" x 8.50", gave a 
diameter-to-width ratio (d/w) of 0.118 and a diameter-to- length 
ratio (d/l) of 0.100. A comparison is made in Section III.C.3 
between the solutions for finite and infinite plates of otherwise 
equal thickness and material constants. 



34 



A fabric G/Ep prepreg material was chosen because it has been 
somewhat less studied than tape and because it is finding 
increased use in airframe construction. The cured fabric laminate 
has slightly less in-plane stiffness and strength per unit thick- 
ness than uniaxial tape made from identical fibers. This is due 
to the nature of the weave, where the fibers (or tows) are cured 
with "crimps" rather than straight. Fabric has, however, demon- 
strated significant advantages over tape in its damage tolerance 
[Ref. 28] and ease of manufacture [Ref. 29]. 

Graphite/ epoxy unidirectional tape can be most effectively 
applied in flat or slightly curved structures such as wings and 
access panels. Fabric, on the other hand, lends itself to appli- 
cations requiring high curvature or complex shapes. Tape cannot 
be used in small inside or outside radius applications without 
fiber separation, inducing matrix- rich/ fiber-poor areas and 
suffering severe loss of strength. 

HMF330C/34 faerie G/Ep manufactured by Fiberite Inc. was 
chosen because it is a high modulus fabric, using Thornel T300 
graphite fibers, found in many aerospace applications. It is an 
eight harness satin (8HS) weave cloth which minimizes the number 
of fiber crimps while maintaining many of the desirable character- 
istics of cloth. Figure 2.2 illustrates some details of its weave. 

1. Panel Reinforcement Conf igura tion 

Reinforcement of the panel cutout was of three general 
types: round, square, and strip. The round and square were 
concentric with the hole, the "stiffening" strips were centered 
0.750 inen away from the hole edge, parallel to the applied load. 
Table I lists the panel designations, reinforcement geometries and 
amounts; Figure 2.3 shows representative configurations. 
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WARP END 




FLOAT 




Figure 2.2 8 Harness Satin (8HS ) Weave Cloth. 
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The basic panel was a quasi-isotropic eight ply (nine 
separate layers including the core) G/Ep panel. For more simple 
comparison, the amount of reinforcement was normalized by the 
amount of G/Ep removed from the 1.00 inch diameter cutout in the 
facesheets of the unreinforced panel. The relative volume of the 
reinforcement ply(s) was determined from this volume (0.088 in ) 
of G/Ep. The round and strip reinforced panels had 5 increments 
of 100% of the removed reinforcement volume and the square rein- 
forcement had increments of 100, 300 and 500%. The 200% and 400% 
reinforcements were each two plies thick. 




Figure 2 . 3 Panel Reinforcement Configurations . 

The panel designation was devised to be somewhat descrip- 
tive of the test specimen. The first letter, ? or R, refers to 
either a £lain (unreinforced) or reinforced configuration, respec- 
tively. The second letter indicates the type of reinforcement: 
none (O), round (R), square (S) or strip (H); X indicates no hole 
was present. The first numeral represents the normalized percent 
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of reinforcement, 1 to 5 for 100% to 500% (0 indicates no rein- 
forcement). The second numeral is the number of reinforcing plies 
on each facesheet. For example, RH42 is a reinforced panel with 
four times the removed hole volume (0.352 in^) arranged in a strip 
configuration, 2 plies thick on each facesheet. 



TABLE I 

TEST SPECIMEN MATRIX 



Panel 

Designation 




Reinforcement 


Type 


Normalized 
Volume ( % ) 


Ply(s) per 
Facesheet 


PX00 


No cutout or reinforcement 


PO00 


None 


0 


0 


RR11 


Round 


100 


1 


RR22 


Round 


200 


2 


RR31 


Round 


300 


1 


RR42 


Round 


400 


2 


RR51 


Round 


500 


1 


RS11 


Square 


100 


1 


RS31 


Square 


300 


1 


RS51 


Square 


500 


1 


RH11 


Strip 


100 


1 


RH22 


Strip 


200 


2 


RH31 


Strip 


300 


1 


RH42 


Strip 


400 


2 


RH51 


Strip 


500 


1 



Figure 2.4 shows a typical laminate cross-section from the 
midplane. Each panel was symmetric about all three axes. 
Exceptional care was required and taken during the manufacturing 
process to ensure that the reinforcement plies were placed 
directly opposite each other on the opposing facesneets. When 
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measured, no reinforcement was more than 0.05" off center; the 
average was less than 0.02". 

D. SOME CHARACTERISTICS OF NOTCHED GRAPHITE/ EPOXY PLATES 

The characteristics of composite materials differ radically 
from those of the metals they replace. As previously noted, 
composite fibers, particularly G/Ep, are by nature very brittle. 
Tensile failures in composite plates with cutouts are, almost 
without exception, load dependent. [Refs. 20-22] 




Compressive failure, the type dealt with here, is more 
dependent on the type and thickness of the laminate, the use of 
honeycomb to overcome the tendency to buckle, the size of cutouts 
and the presence of imperfections. The compressive failure modes 
tend to be complex, composed of one or more types of failure: 
stability, ply delamination, matrix cracking, etc. Stability 
failure is principally the buckling of either fibers within the 
matrix (micro-mechanical) or the structure itself (macro- 
mechanical). These test specimens and the frame were designed to 
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preclude macro-mechanical buckling (in the Euler column mode) 
since little would be learned about the reinforcement effects and 
this type of failure has been well documented beginning with Levy, 
Woolley and Kroll [Ref. 5]. 

1. Stress Concentration due to Notch Effects 

It is well known that notches and cutouts in plates act as 
stress risers. For circular holes in plane elastic isotropic 
infinite plates under uniaxial tension or compression, the stress 
at the hole edge 90° to the applied load will be exactly three 
times the far-field stress. The distribution of stress around the 
hole edge and the stress field around it can be predicted using 
Airy's stress function. Dally and Riley [Ref. 30: pp. 67-83] give 
a clear and concise derivation of the stress field equations which 
will not be repeated here. 

2. Orthotropic Effects on Stress Distribution 

When an orthotropic plate with a stress riser is loaded, 
the SCF depends on the degree of orthotropy, that is, how much the 
elastic modulii change with radial direction. This is sometimes 
referred to, not always correctly, as the ratio of E 1 /E t . The 
subscripts "1" and "t" refer to the effective lateral and trans- 
verse modulii where the lateral direction is parallel to the 
applied load and transverse is 90° to it. In the coordinate 
system used in this report, the load is applied parallel to the Y 
axis and the ratio is expressed as: Ey/E x * Note that a ratio of 

1.0 does not ensure isotropy; it must be accompanied by the appro- 
priate shear modulus (G„.,) and Poisson's ratio ( £>,)• For a 

A y A y 

circular hole in an infinite-width plane orthotropic plate, the 
stress concentration K M on the cutout edge 90° to the applied 
load was given by Nuismer and Whitney [Ref. 27: Eqn. 3] as: 
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K 



1 + 



( 2 . 1 ) 



/ 2< JE l/E t - t’lt) + Ei/Git ' 

In an idealized infinite laminated plate, this equation 
must be equally valid in both tension and compression. This 
stress concentration factor (K ) may be considered a far better 
indication of the orthotropy of a material than the ratio E^/E t . 

The distribution of stress in the Y direction along the X 
axis ((Ty(x,0)) due to an applied (far-field) normal stress (o" n ) 
may be approximated using the following equation: 

(7 y (x,0) = [( a n /2)(2+b 2 +3b 4 -(K -3 ) (5b 5 -7b 8 ) ], (2.2a) 

where : 

b = a/ (x-d/2) and x > d/2. (2.2b) 

The variable "d" is the diameter of the circular cutout and "x" is 
a location along the X axis (y = 0.0") when the coordinate system 
is concentric with the hole. This relationship is a quite 
accurate polynomial approximation developed by Konish and Whitney 
[Ref. 31]. 

3. Effects of Finite Plate Width and Hole Size on SCF 

Compared to infinite plate width under uniaxial stress, 
finite plate width acts to increase the SCF. This fact becomes 
obvious in plates with a high d/w ratio. The applied stress must 
be carried by a greatly reduced net cross-section. The increase 
in SCF is due more to the net section effect than the presence of 
the cutout. Peterson [Ref. 32: pp. 110-111] gives the following 
equation to approximate the SCF at the edge of an unreinforced 
circular cutout in a finite-width isotropic plate: 

K f = [2 + (l-(d/w) 3 )]/[l-(d/w)] (2.3a) 
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This can be extended to an orthotropic plate where does not 
equal exactly three using: 

K f = ( K ®/ 3) 2+[l-(d/w) 3 ] /[l-(d/w) ] (2.3b) 

The test specimens used in this report had a d/w ratio of 
0.118; Kf was then calculated to be 3.045 for the unreinforced, 
quasi-isotropic panel PO00. At TT n = -10.0 ksi this would 
theoretically make the maximum stress -30,450 psi at the Q = 0° 
position on the cutout (90° to the applied load) compared with 
30,000 psi predicted for an infinitely wide plate. This is an 
increase of 1.5%. Thus panel width has little more than a 
negligible effect on the SCF of the test specimens in this report. 

Further data that relate a plate's dimensions to its SCF 
are given by Hong and Crews [Ref. 25: pp. 8-10]. They calculated 
stress concentration factors in finite-width orthotropic plates 
under uniaxial loads using finite element analysis. They used a 
different definition of SCF, one based on the net cross-sectional 
area stress concentration (K ne f). This report uses the SCF based 
on far-field stress or the gross SCF ( K g ross )* The two are 
related by the equation: 

V°ss - K net /[l-(d/»)]. (2.4) 

To make valid comparisons with SCF data presented in this report 
selected results of Hong and Crews' analysis, converted from K nat 
to K g ross , are listed in Table II. 

Their results show that quasi-isotropic layups 
(C0,±45,90] s ) give results very close to the theoretical isotropic 
values. Greater orthotropy in the load direction results in a 
correspondingly greater SCF. It is interesting to note that the 
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ratio of length to width (1/w) has an increasing effect on the SCF 
as the ratio d/w increases. 

Nuismer and Whitney [Ref. 27: p. 118] point out the effect 
of absolute hole size on panel failure: ". . . attention was 
called to a phenomenon that since became known as the 'hole size 
effect, 1 that is, for tension specimens containing various sized 
circular cutouts, larger holes cause greater strength reductions 
than do smaller holes." They state that the classical stress 
concentration approach does not explain such an effect and they go 
on to propose that while the stress concentration factor is the 
same, the distribution and gradient near the hole is different. 
Figure 2.5 reproduced from Ref. 27 illustrates this point. 

4. Failure Stress Criteria 

As previously noted, the SCF does not explain the "hole 
size effect" on failure. Nuismer and Whitney rejected linear 
elastic fracture mechanics (LEFM) to explain the inverse relation- 
ship between hole size and strength. They noted that while all 
circular holes in infinite width plates should have the same 
theoretical SCF, the distribution in fact changes with hole 
radius. The smaller the hole the more concentrated the stress 
near the edge appears [Ref. 27: p. 118]. Nuismer and Whitney 
proposed that when the notched stress (CT^) reached an average 
value of (J u>ult , the unnotched ultimate stress over some 
characteristic distance ao, that the panel's ultimate strength had 
been reached and failure resulted. This characteristic distance 
ao must be arrived at by testing a statistically significant 
number of panels. This distance ao is defined: 



TABLE II 



STRESS CONCENTRATION FACTORS-K gross (HONG & CREWS) 



Layup 




El /E t 


K L/w 


0. 


Diameter-tO' 
.05 0.10 


-Width 

0.33 


Ratio: 

0.50 


d/w 

0.67 


[0,±45, 


90] s 


1.00 


3.00 


1 




3.00 


2.74 


3.33 


4.02 


5.76 


[0 , +45 , 


90] s 


1.00 


3.00 


2 




3.01 


3.03 


3.49 


4.36 


6.36 


[90] 




0.07 


2.48 


1 




2.48 


2.51 


2.97 


3.78 


5.61 


[90] 




0.07 


2.48 


2 




2.48 


2.51 


2.97 


3.78 


5.88 


[±45] s 


★ 


1.00 


2.06 


1 




2.88 


2.93 


3.38 


3.84 


5.16 


[±45] s 


★ 


1.00 


2.06 


2 




2.88 


2.92 


3.36 


3.80 


5.28 


[0 , 90 ] s 


★ 


1.00 


3.78 


1 




4.78 


4.69 


5.61 


5.08 


6.93 


[0,90] s 


★ 


1.00 


3.78 


2 




4.82 


4.84 


5.22 


6.06 


8.16 


1 — 1 

O 
* j 




13.49 


5.43 


1 




6.36 


6.07 


5.24 


5.82 


8.01 


i — i 

0 

1 1 




13.49 


5.43 


2 




6.44 


6.44 


6.54 


7.30 


9.54 




* Indicates an E,/E t 


ratio 


of 1 . 


00 , but 


not a 







quasi-isotropic laminate. 




(Reproduced from Reference 27) 

Figure 2.5 Hole Size Effect on Normal Stress Distribution 
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The ratio of notched to unnotched ultimate strength (O n /O u ) u i~ 
for infinite plates is: 



° n /0 



U 



where : 



2(1- £ ) / C 2— £ 2 - £ 4 +(k q 3 -3)(| 6 - £ 8 )] 
£ = d/(d/2 + a«.) and x > ao. 



(2.6a) 



(2.6b) 



Nuismer and Labor [Ref. 33: p. 55] determined that for 
AS/3501-5 G/Ep (tape) in compression this characteristic length 
was 6.2 mm (0.24"). They also note that the characteristic length 
for tape in tension was only 2.3 mm (0.091"). Test data provided 
by LMSC indicates that for HMF330C/34 fabric G/Ep this character- 
istic length is close to 7.3 mm (0.33") in compression. 

5. Effect of Poisson and Interlaminar Stresses on Failure 
Isotropic materials may be modeled using classical plate 
theory neglecting out-of-plane stresses (±z in this coordinate sys- 
tem). Orthotropic materials, however, develop complex inter- 
laminar stress fields near the edge of a cutout. The subject has 
received much theoretical attention [Refs. 34 through 36]. Tang 
[Ref. 34: p. 1631] states that ". . . radial and shear stresses of 
each layer along the contour of the hole are in general not zero 
because there exists a three-dimensional state of stress at the 
free edge of each layer which the plane stress solution cannot 
predict." Greszczuk [Ref. 9: p. 372] pointed out that "In ortho- 
tropic and anisotropic plates containing openings, the failure 
will take place not as a result of stress concentration, but 
rather as a result of interaction of various stress components." 

Under uniaxial compressive loading the laminate will have 
a Poisson expansion induced out-of-plane tensile stress (<X Z ) 
which is highest at the hole's edge at point of the greatest 
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stress concentration. This stress is added to any local stress 
due to machining and imperfections and combined tend to hasten 
delamination and the ultimate failure. In the experimental 
results reported here it was not possible to effectively quantify 
the effect on failure of this out-of-plane stress. 
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III. COMPUTATIONAL ANALYSIS 



Before an experimental program could be developed, it was 
necessary to understand and be able to analyze the strain field 
resulting from a cutout in a representative panel and to be able 
to predict the reaction of test specimens to an applied com- 
pressive load. Three analysis methods were used: classical 
laminate theory, the linear elastic stress function and linear 
elastic finite element analysis (LEFEA). 

Laminate analysis provides the basic stress-strain relation at 
a point, once the material properties of each constituent ply are 
specified. The stress function was used to predict the theoreti- 
cal stress-strain fields in an infinite unreinforced orthotropic 
elastic plate with a circular cutout. These two can be solved in 
closed-form and require relatively little computation time using 
modern computers. The finite element method allows detailed 
analysis of reinforced finite-width reinforced panels, but 
requires a significant allocation of computer resources for an 
accurate representation of the strain field. 

There are several coordinate systems and notations in general 
use in laminate analysis. The following section presents the 
method used in this report, explicitly defines the notation and 
gives justification for some of the assumptions that were made. 

A. LAMINATE THEORY AND ANALYSIS 

Laminate theory seeks to predict the properties of a multi- 
directional composite laminate based on the properties and 
orientation of its constituent lamina. Individual laminae are 
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usually either unidirectional (tape) or woven (cloth) fibers 
embedded in a polymer matrix (generally a thermoset resin whose 
molecules are linked in three dimensions and which exhibits 
elastic properties in normal use) and tend to have strongly 
directional properties. Tne theory assumes that the state of 
stress is plane, displacements are small compared to laminate 
thickness and that strain is much smaller than unity. 

Pipes [Ref. 37: pp. 4-1, 5-1] presented the micro- and macro- 
mechanical models tnat are the basis of the theory. An anisotro- 
pic material's elastic response at a point to applied stresses may 
be defined using generalized Hooke's Law. The constitutive rela- 
tion is Equation 3.1, where CJj_j and 6 ^ are the components of 
the stress and strain tensors and C— ^ is the tensor of elastic 
constants. Using this most general of equations there are 3^ or 
81 material constants. 

3 3 

^ij “ ^ > C ijkl ^kl (3.1) 

k=l 1=1 

This equation may be greatly simplified using the symmetry of 
stress and strain and the requirement that the strain energy 
density function be positive definite [Love, Ref. 38: pp. 97-111 
and Feynman, Ref. 39: v. 2, ch. 31-7] reducing the independent 
elastic constants from 81 to 21. Symmetry reduces both <T and 
£- K ]_ from nine to six different values. Feynman explains that the 
elastic response of a crystal with no symmetry in the three axes 
can be completely defined using 21 independent coefficients. The 
notation can be contracted using the following convention where 
the index: i = 1,2,3: 
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(3.2) 



li 








The constitutive relation can now be expressed as: 



-E 

j=i 



1 3 



(3.3) 



This is a sixth-order symmetric matrix (where C^j = Cjj_). 

Idealized thin laminate theory neglects stress and strain in 
the ±z direction; the equations are reduced to plane strain and 
stress, further contracting the elastic constant tensor to a 
third-order symmetric matrix. In orthotropic systems (axes at 
right angles to each other) the "1" direction is the principal 
fiber direction (or the direction with the highest elastic modu- 
lus), "2" is 90° to it and "6" is the shear in the 1-2 plane. 



°i = Q ij £ j’ 

j— 1, 2,6 



(3.4) 



The matrix [Q] is termed the reduced laminate stiffness matrix, is 
symmetric and is related to [C] by: 

c . c 

Qj_j = C ± j — (i,j = 1,2,6) (3.5) 

C 33 

The matrix Cq] may be expressed explicitly in terms of modulii 
and Poisson's ratios: 



Qll = 



1 “ ^12 ^21 



Ql2 “ 



^12 E 2 



\ - V V 

1 12 21 



(3.6) 
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(3.6) 



S 2 


II 

i — 1 

(N 

O 


Ql 2 = 0 


22 " 1 - "l2 V 21 


Q 16 “ 


Q 61 = 0 


66 = G 12 


Q 26 = 


Q 62 = 0 



To determine the ply's elastic response defined in the 
laminate coordinate system (X-Y), both the stress and strain 
vectors must be rotated an angle Q about the "3" axis (Note: the 
"3" and Z axes are colocated): 



^ T aijCJj (i = x,y,xy) (3.7) 

j= 1/2,6 

C-i = T ei j £j (i = x, y , xy ) (3.8) 

j= 1,2,6 

The transform matrix [T a ] (for the case of stress) is derived from 
the trigonometric relations: 



m 2 


n 2 


2mn 


m = cos 




n 2 


m 2 


-2rnn 




(3.9) 






0 7 


n = sin 




-mn 


mn 


m -n 2 







Recall tnat engineering shear strain (T^) differs from 
tensorial shear strain ( £ ^) by a factor of 2: 6- = 2T 6- The 

strain transform matrix elements T^j^ and T^g become ran and T*^ 
and Tf ^2 become 2mn. Using matrix algebra, the now transformed 
reduced laminate stiffness matrix [Q'J can be expressed in matrix 
form as : 
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[Q*] = [T a ] [Q] [T e ] 



(3.10) 



-1 



A laminate is built from the stacking of a number of these rotated 
plies. The designer may easily tailor the laminate using various 
ply thicknesses and orientations. 

The integration of each ply's [Q'3 matrix through the laminate 
thickness (h) gives the normalized inplane stiffness matrix: 



[A] = 




(3.11) 



Stress and moment resultants are defined by integrating stress 
through the laminate thickness: 




(3.12) 




(3.13) 



The stress resultant vector is related to the strain 

vector j £ | by the laminate inplane stiffness matrix [A] in 
equation 3.14: 



N x 




' 1 

' 1 

< 

1 


A 12 


A lb' 


i €X 


Ny 


^ = 


A 21 


A 22 


A 26 


£y 






1 

> 

(T\ 

I— 1 


A 62 


A 66. 


( ^ xy 



(3.14) 



The laminate, while thin, demonstrates resistance to bending 
governed by the ply stiffness and the square of the distance from 
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